function import_xyz(obj,varargin)
%Select atoms of a given element


fid = fopen(char(varargin{1}),'r');

natom = fscanf(fid,'%i');


for ka=1:natom
    
    kind={fscanf(fid,'%s',1)};
    position(1)=fscanf(fid,'%f',1);
    position(2)=fscanf(fid,'%f',1);
    position(3)=fscanf(fid,'%f\n',1);
    obj.add_atom(position,kind);

end

fclose(fid);    
end
